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Oh| Abstract: Recently, an exact description of instanton corrections to the moduli spaces of 

r-| \ 4d = 2 supersymmetric gauge theories compactified on a circle and Calabi-Yau compacti- 
fications of Type II superstring theories was found. The equations determining the instanton 

^ ^ contributions turn out to have the form of Thermodynamic Bethe Ansatz. We explore further 

\ this relation and, in particular, we identify the contact potential of quaternionic string moduli 

^ \ space with the free energy of the integrable system and the Kahler potential of the gauge theory 

\ moduli space with the Yang- Yang functional. We also show that the corresponding S-matrix 

g \ satisfies all usual constraints of 2d integrable models, including crossing and bootstrap, and 

\ derive the associated Y-system. Surprisingly, in the simplest case the Y-system is described by 



the MacMahon function relevant for crystal melting and topological strings. 
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1. Introduction 

The last years have marked a remarkable manifestation of integrability and in particular of 
(Thermodynamic) Bethe Ansatz in description of dynamics of gauge and string theories. A 
prominent example is the calculation of the spectrum of conformal dimensions in = 4 super 
Yang-Mills theory [|^ and the exact solution of its AdS / CFT dual given by a non-linear a- model 
^ . Nowadays this is a wide area of research where the integrability in the form of spin chain 
models and Bethe Ansatz plays the crucial role. However, during the last year several new 
connections to the integrable world have emerged. 

First, Nekrasov and Shatashvili discovered 0, ^ that N = 2 supersymmetric gauge theories 
are in one to one correspondence with integrable Hamiltonian systems in such a way that the 
supersymmetric vacua of gauge theories are mapped to Bethe states of the integrable models. 
Thus, their identification heavily relies on Bethe Ansatz and its appearance in the description of 
BPS vacua. Another interesting development has been done by Alday, Gaiotto and Tachikawa 
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who related the instanton partition functions of certain N = 2 SCFTs to conformal blocks 
of Liouville theory. As was recently demonstrated in the latter duality is actually closely 
related to the former one. 

On top of that, it was realized in ^ (see also [|, |lOl) that the classical problem of deter- 
mining the minimal area surface in AdS^ ending on a null polygonal contour, which provides 
the strong coupling limit of gluon scattering amplitudes in iV = 4 super Yang-Mills [^, has 
a solution in terms of Thermodynamic Bethe Ansatz (TBA) |jl2[. The TEA free energy then 
gives the area encoding the scattering amplitudes. 

Here we consider another example of the unexpected interplay between gauge/string theory 
and integrability whose physical explanation is still lacking. Namely, it was noticed in |TB[ (see 
also jl^) that equations describing the exact moduli space of 4d = 2 supersymmetric gauge 
theory compactified on a circle coincide with the equations of TBA. Moreover, exactly the 
same equations appear also in the twistor description of D-instanton corrected moduli space 



of Type II string theory compactified on a Calabi-Yau (CY) threefold ]T5|, |T6|. Due to this 



one may hope that the elaborated machinery of integrable systems will give new insights for 
these research areas as it did in the above mentioned situations. However, for this one needs 
to extend the correspondence beyond just the level of equations to some interesting physical 
quantities as well as to understand better the integrable structure behind this TBA. 

In this paper we try to fulfil this goal by pushing forward the relation between TBA and the 
instanton corrections to the moduli spaces. In particular, we identify the TBA free energy with 
the instanton contribution to the so called contact potential governing the quaternion-Kahler 
geometry of the string moduli space. In the gauge theory context the analogous quantity is 
the Kahler potential. We show that its most non-trivial part turns out to be encoded by the 
Yang- Yang functional [17] of the associated integrable system. 

To uncover the integrable structure, we study the S'-matrix which can be found from 
the TBA equations. Although it is not unitary |jl3[, we show that it satisfies all standard 
constraints usually imposed on factorizable S-matrices of integrable systems. They include 
Lorentz invariance, crossing symmetry, Yang-Baxter equation and bootstrap identity. Finally, 
we derive the Y-system following from our TBA equations. 

The paper is organized as follows. In the next section we review the twistor description 
of hyperkahler and quaternion-Kahler spaces which is crucial for introducing general instanton 
corrections to the moduli spaces. These instanton corrections are further described in section 
2.2 where we present the equations determining the non-perturbative geometry of the moduli 
spaces and the results for the contact and Kahler potentials. To our knowledge, the last quantity 
given in ( ^.31| ) was not known before. The equations ( |2.25| ) represent our main starting point 
and the reader not interested in their origin can pass directly to them. In section ^ we establish 
the relation of (|2.25| ) to TBA and express the contact and Kahler potentials in terms of the 
free energy and Yang- Yang functional, respectively. After that in section | we investigate the 
S'-matrix associated to our problem, derive the Y-system and comment on the conformal limit 
of our TBA. In section ^ we conclude with a discussion and some open questions. In appendix 
^ we provide some details about evaluation and symplectic invariance of the Kahler potential. 
Finally, appendix O is devoted to a particular case of rigid CY where we observe the appearance 



of the MacMahon function 18 
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2. Non-perturbative moduli spaces in gauge and string theories 



2.1 Twistor description of HK and QK spaces 

The moduli spaces of compactified gauge and string theories with N = 2 supersymmetry 
are examples of manifolds with a quaternionic structure. The supersymmetry restricts the 
moduli space to be hyperkahler (HK) in gauge theory [|I^ and its hypermultiplet sector to 



be quaternion-Kahler (QK) in string theory ^0 . In both cases the geometry of these moduli 
spaces gets perturbative and non-perturbative contributions and, of course, it is extremely 
difficult to describe them explicitly as corrections to the metric. It would be much more 
convenient to have at our disposal some kind of prepotential, which can encode all quantum 
corrections in a systematic way, similarly as does the holomorphic prepotential F{X) for the 
vector multiplet moduli space. Thus, we are confronted with the problem of parametrizing 
quaternionic geometries. 

This problem is solved by considering the twistor space Z_\4 of a quaternionic manifold 
Ai, which is a CP^ bundle over Ai and provides its very efficient description. In particular, it 
allows to encode all geometric information in a set of holomorphic functions playing the role 
of prepotentials we were asking for. Here we briefly review the corresponding construction and 
refer to ^ for more details. 

The twistor description is so powerful because of the existence of a certain holomorphic 
structure on Z_m. In the HK and QK cases, these structures are represented locally by holomor- 
phic forms, a 2- form fl and a 1-form X, respectively. The former deflnes a symplectic structure 
on each fiber of the bundle n : Zm — )■ CP^, and the latter produces a contact structure on 
Zj^4. As usual, locally it is always possible to choose Darboux coordinates where these forms 
take a standard form. In other words, the twistor space can be covered by a set of open patches 
Ui such that in each patch one has 

HK: ^]['']=d/il;Ud^/f], QK : X^"^ = da^^^ + ^^^d^f , (2.1) 

where A = 0, 1, ... ,6? — 1 with Ad = dim^ A^, and {u^-^, jJ^ , () (or (^pj; ' form a set of 
holomorphic Darboux coordinates in the patch Ui of the twistor space of HK (respectively, QK) 
manifold. Here C is a complex coordinate on <CP^ (also in the QK case) and we will use the 
notation to parameterize the base quaternionic manifold Ai. 

In addition to this holomorphic structure, the twistor space carries a real structure defined 
in term of the antipodal map r acting on CP^ as t : ( ^ To express the compatibility 

of the two structures, we assume that the antipodal map sends Ui to Wj. Then the holomorphic 
forms should satisfy the reality constraint 



r 



(^]H) = ^]^^ r(A'H) = A'l^ (2.2) 



so that the Darboux coordinates may be chosen to satisfy similar relations, all with sign plus, 
under the combined action of the complex conjugation and the antipodal map. 

Something non-trivial appears when one considers the overlap of two patches Ui (lUj. In 
this region the holomorphic forms defined in different patches must be related as 

fil^l = 4 mod dCt^ A'W = A'bl, (2.3) 
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where fij are (fixed and simple) transition functions of the 0{1) bundle on CP^, whereas fij 
are fixed later. This implies that the local Darboux coordinates are related by symplectic 
and contact transformations, respectively. Such transformations are generated by holomorphic 
functions on Z_m, which we call transition functions H^^^^. It is these functions that play the 
role of the prepotentials. They establish relations between Darboux coordinates in different 
patches, which together with the reality constraints and suitable boundary conditions allow to 
find these coordinates as functions of ( and x^. These solutions, called twistor lines, contain 
all geometric information and are the starting point to find the metrics on Ai and its twistor 
space. 

In the following we restrict our attention only to those transition functions which in the 
QK case are independent on the coordinates a'*'. Physically, for the hypermultiplet moduli 
space of compactified Type II string theory this means that we ignore the contributions due 
to NS5-brane instantons [^. Such restriction brings great simplifications. In particular, the 



coefficients fij left undetermined above are actually given by 

i;5 = i-9„,ii/[^^-i (2.4) 

and thus reduce to 1 in this case. What is important for us is that given this restriction the 
equations determining the twistor lines in HK and QK cases become equivalent. Indeed, in the 
QK case for a-independent //^'■'l they look as 

^b] = #1 - = + d^.H^^^ (2.5) 

whereas in the HK case they are the same provided one replaces C,^-^ by r]^^ = {iC)~^ flii'^-^ and 

by iJj^ . Here the index + refers to some fixed patch W+ which we choose to be the one around 
the north pole C = 0- The coincidence of the equations allows us to consider the HK and QK 
cases simultaneously. The notations mostly used in the following correspond to the QK case. 

The gluing conditions ( p.5|) for the twistor lines can be rewritten in a form more convenient 
for a perturbative treatment as the following integral equations ||16|| ^ 

where ( E Ui, Cj is the contour surrounding Uj in the counterclockwise direction, whereas 
complex and real A^,B\ are free parameters playing the role of coordinates on A^.^ The 
sum in ( p.6[) goes over all patches including those which do not intersect with Wj. In that case 



the transition functions are defined by the cocycle condition and by analytic continuation |£l 



"'^ Comparing to |13, |5|, we changed some normalizations to avoid some numerical factors and to make the 
symplectic invariance more explicit. In particular, /ha, ^a, a and ijl*-'! are renormalized by factor — 2i, f2, e'^ 
by 2, and i^^ by i. 

^Actually in the QK case the phase of F° can be absorbed into ( so that it becomes real. The missing 
coordinate is provided by the real part of the constant coefficient Ba in the expansion of a^~^\ which we do not 
write here explicitly. 
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Once these integral equations are solved, there is a straightforward procedure to extract 
the metric p2| , |T^. On this way the prominent role is played by the Kahler potential Kj^^ in 



the HK case and the so called contact potential in the QK case, which provides a Kahler 
potential for the metric on the twistor space Zj^^ 



i^| = logi^+ Re$H. 



(2.7) 



The contact potential appears as a certain coefficient in the expansion of the contact one-form 



A'H=2-^(dC+p+-ip3C+P-C') 



(2.8) 



where p is the SU{2) part of the Levi-Civita connection on Ai, and generically it is a func- 
tion <l>[j](x^,C) holomorphic on the CP^ fiber, but defined only locally what is refiected by 
the patch index. However, in the case under consideration when the transition functions are 
a- independent, it becomes real, globally defined and independent on (, = 0(a;^). The 
remaining function on the base manifold is given by 



(2.9) 



where is a constant called anomalous dimension, which encodes a boundary condition for 
the twistor line aW. 

In the HK case the Kahler potential can also be expressed as an integral of transition 
functions. In our notations the representation found in [|16| can be summarized as 



K 



M 



4^ 



(2.10) 



By simple manipulations using the integral equations (|2.6|) , it can be rewritten in the following 
form 

^M = iT. £ J [^"' - 4AtH''' + {(-'y' - cy"-) ^fH^'^'] ■ (2.11) 

This shows that in the physically important case of transition functions homogeneous of degree 
one, relevant when the HK space is the hyperkahler cone of a QK space [H, 22|, the Kahler 
potential essentially coincides with the contact potential (p.9|). 

We close this subsection by noting that although all integration contours Cj appearing 
in the formulae above are closed since they surround open patches, under some conditions 
it is possible to generalize these results to include also open contours. Typically, such open 
contours can be viewed as a leftover of some discontinuity cuts in transition functions of a 
"more fundamental" construction using only open patches and closed contours. The former 
arises from the latter when one shrinks the contour surrounding the cut so that its contribution 
reduces to an integral of the discontinuity of the initial integrand. Such description in terms 
of open contours turns out to be relevant in the discussion of the instanton corrected moduli 
spaces in the next subsection. 
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2.2 Instanton contributions to the moduli spaces 



The low energy dynamics of = 2 supersymmetric gauge and string theories is completely 
determined by the geometry of their moduli spaces. For a 4d iV = 2 gauge theory with the 
gauge group G of rank rf, which is compactified on a circle of radius R, the moduli space 
is parameterized by complex scalars from the vector multiplet and by Wilson lines of the 
gauge potential around the circle, which have "electric" and "magnetic" components Ca 
and are all periodic.^ It is a hyperkahler manifold whose perturbative metric follows from the 
simple 3d truncation of the 4d vector multiplet Lagrangian and thus is completely defined by 
the holomorphic prepotential F{z). However, it gets instanton contributions from the massive 
spectrum due to BPS particles going around the compactification circle. For large R these 
contributions are exponentially suppressed since for a particle of charge 7 = {qA,p^) they are 
weighted by c'^'^-^I^tI where 

Z^{z) = qAz''-p''FA{z) (2.12) 

is the central charge function giving the mass of the BPS particle. 

The low energy physics of Type IIA string theory compactified on a CY threefold X is 
a bit different. Its complete moduli space is factorized to the moduli spaces of vector and 
hypermultiplets. The former is tree level exact so that our interest is concentrated on the 
hypermultiplet sector. It comprises d = h2^i{X) + 1 hypermultiplets, which include the complex 
structure moduli = J^^ Fa = J^^ Q, the RR scalars Ca representing the RR three-form 
integrated along a symplectic basis (7^, 7a) of A and B cycles in H^^X, Z), the four- dimensional 
dilaton e*^ = l/5'^4) and the Neveu-Schwarz (NS) axion a, dual to the NS two-form B in four 
dimensions. Whereas X^ provide a set of homogeneous coordinates for complex structure 
deformations, they may be traded for the inhomogeneous coordinates z"" = X"/X°. As we 
mentioned above, this moduli space is a quaternionic-Kahler manifold. Its tree level metric 
is determined by the holomorphic prepotential F{z) as in gauge theory and receives one-loop 
P3| , PU| and instanton corrections [^. The latter arise either from D2-branes wrapping non- 
trivial 3-dimensional cycles of X or from the NS5-brane wrapping the whole Calabi-Yau. In 
this work we ignore the second type of instantons so that we remain only with membrane 
contributions. Since cycles in H^i^X, Z) are parameterized by the symplectic vector 7 = {qA,P^) 
and the weight of the corresponding instanton is determined by the same function (|2l2|) (with 



z° = 1), one may expect that the D-instanton corrected hypermultiplet moduli space and the 
exact gauge theory moduli space from above have a similar description. 

This indeed turns out to be the case. Such description is provided by the twistor formalism 
from the previous subsection, which as we have seen works in the same way for HK and QK 
spaces. Below we review the corresponding construction for the gauge and string moduli spaces 



ISj 0. It essentially amounts to provide a covering of CP^ and a set of transition functions. 
First, we present it for the moduli spaces at the perturbative level and then show how the 
instantons are included. 



■^We use different notations than those of The fields are denoted by the same letters which are used to 
parametrize the moduli space of Type IIA string theory since in all equations they appear in exactly the same 
way. Besides, the Wilson lines are normalized to have period 1. 
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To construct the perturbative moduli space, let us cover the Riemann sphere by the follow- 
ing three patches: the first patch W+ surrounds the north pole, the second patch U_ surrounds 
the south pole and the rest is covered by Uq (see Fig. |l|a). There are two non-trivial transition 
functions associated with such covering. 

ifM=F(e[+]), //[-°]=F(eH)- (2.13) 

As usual, in the gauge theory case one should replace ^[j] by //[jj.^ On the other hand, in the 
string theory case one should add information about the anomalous dimension (see (^^) ). 
It turns out to be determined by the Euler number of CY, Cq = xx/(967r), and incorporates 
the one- loop correction p2| . Then the equations ( p.6|) provide an explicit representation for the 



twistor lines. Identifying properly the abstract coordinates Z'^, A^, Bf^ of the twistor approach 
with the physical fields, in the patch Uq (in string notations) one finds 



efo, = efpe.)-CV7^(^V-c^^), 

if = = Ca + 7^ {(-'Fj^iz) - c Fa(^)) , 

where TZ = R/2 va. the gauge theory case and can be expressed through the dilaton, using 



with K[z,z) = — 2Im(^^FA), in string theory. Remarkably, the contact potential ( |2.9| ) turns 
out to coincide with the dilaton, which is in turn related to the four- dimensional string coupling 
(7(4). On the other hand, TZ is inversely proportional to the ten-dimensional string coupling (7(10) • 
These identifications will survive the instanton corrections. Altogether the above results are 



sufficient to extract the metric in a straightforward way |22 . 

Let us include the instanton contributions from a set of 2N BPS particles (or D2-branes) 
with charges {7a}a=i- The number of particles is taken even because every particle is accom- 
panied by its antiparticle with the opposite charge. Eventually we are interested in the limit 
iV — )■ 00. We assume that the charges are ordered in accordance with decreasing of the phase of 



(|2.12|) and that these phases for the charges 7^ which are not mutually local, i.e. for those 



which have non-vanishing symplectic invariant scalar product 

(71,72) = gi,AP^-g2,AP^ (2.18) 



^In fact, in the gauge theory case the holomorphic prepotential contains a logarithmic term and thus it is 
only a (/Masi-homogeneous function of second degree 

z^Fa^2F+^Qa^z^z'', (2.14) 

where Qah is a matrix constructed from charges of hypermultiplets whose explicit form will not be important in 
the following. On the other hand, the tree level transition functions should not contain logarithmic singularities. 



Due to this, the functions (2.13) should be actually replaced by 

^[+0] ^ ^-^FiC^^^), = C'F{C\^]). (2.15) 
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Figure 1: Covering of CP and transition functions of the perturbative and instanton corrected 
twistor spaces. 



are all different. For mutually local charges with coinciding phases of the order is not 
important. 

Each charge vector 7 defines a "BPS ray" on CP^ going between the north and south 
poles as 

£, = {C: Z(7)/CGiM"}. (2.19) 

These rays split the patch into IN sectors which we call Ua-, so that the covering of CP^ 
consists now of these connected parts divided by the contours and the usual patches U.± 
(see Fig. |T|b). 

What we have to provide is the transition function through the BPS ray associated with a 
charge 7a. For that purpose we define 



\ah\ 



7c 



^cAil-v'-M- (2.20) 



Then the transition function through the BPS ray is [|T6 



^[aa+l]^^^^j^|[a+l]) ^ - ^ g.,AP^(G;J^ (2.21) 

where (^-^^(S^^) with = si^"' is defined in term of the dilogarithm function as 

G,(S,) = ^Li,(e-^-^). (2.22) 

The coefficient n-y encodes the spectrum of the gauge theory or carries a topological information 



about Calabi-Yau. In the former case it was identified with the second helicity supertrace in |13 



and in general it can be related to the generalized Donaldson-Thomas invariants introduced in 
p7[ . For us two facts are important: are integer numbers and do not depend on the sign of 



the charge, = ri-^. The presence of the second (non-symplectic invariant) term in ( |2.21| ) is 
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related to the fact that the arguments of transition functions should be Darboux coordinates 
from different patches. Therefore, the r.h.s. of (|2.21|) should be expressed through fjf^^^. This 
is achieved by solving 

^r^''=2,.-ga,ApX(S7J (2.23) 

for which, of course, cannot be done explicitly. Nevertheless, this is sufficient to define the 
derivatives of which lead to the following simple symplectic transformations 

efa+i] = efa] + p^g; , i^'^ = + ^..ag; . (2.24) 

To complete the construction we need also to provide if'^^l. They are given by (|2.13| ) plus 
instanton contributions which are known but quite complicated. We refer to [1^ for their 
explicit expressions. 

Note that this construction realizes the option mentioned in the end of the previous sub- 
section. Namely, the transition functions ( p. 21 ) are to be integrated along open contours given 



by the BPS rays. There exists also a version with closed contours |]T5[, which however does not 
bring anything new. 

Another useful comment is that this picture gives an interesting geometric interpretation 
of the wall crossing phenomenon. The lines of marginal stability where it takes place arise at 
the moduli configurations where the phases of the central charge ( |2.12|) of two BPS states align. 



This happens precisely when the two BPS rays (2.19) cross each other. Since in general different 



order of BPS rays leads to inequivalent constructions, this is a clear origin of discontinuities. 
On the other hand, crossing a line of marginal stability, the BPS spectrum, encoded in the 
coefficients n^, also changes discontinuously. A formula recently found by Kontsevich and 



Soibelman ensures that the two effects compensate each other and the resulting metric is 



continuous 13 



Although the presented construction uniquely defines the twistor space and the underlying 
quaternionic manifold, to actually compute the metric one needs to solve the integral equations 
( p.6|) for the twistor lines. In our case, they can be reduced to the following system^ 



5^.(0 = + ^ {C'Z^^ - C^^J + A 5Z (7a, lb) 


C C-C ^ 




(2.25 



where 

G^^^aC^-ZCa. (2.26) 

These equations encode all non-trivialities of the problem. For large TZ (large circle radius or 
small lOd string coupling), they can be analyzed perturbatively and their solution, represented 
by the set of variables S^^(C), contains all orders of the instanton expansion. Note that due 
to the reality conditions on twistor lines and the fact that Ua is the patch associated with the 



^Here we ignore an additional issue of sign in front of the exponent in the logarithm appearing in the gauge 
theory context. It is related to some subtleties in the fermion number of bound states and can be taken into 
account by the so called "quadratic refinement" [|l3| . It affects the solution in a simple way and is not important 
for our discussion. 
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charge ■jg. = —7a, the functions S^^(^) satisfy 



S,.(C) = -5-,„(-n 



(2.27) 



Once the functions S^^(C) are known, the twistor hues easily follow as |]T6 



(.[a] _ ^■(pcrt) 



(2.28) 



where ( and 



^^logfl-e-2-=^(^')V 



C' c-c 

The general formula for the contact potential (|2.9|) gives the following result 



2 ^ ' ^ 967r 167r2 



(2.29) 



E ^> / ^ (C-^^.. - C^>) log (1 - e-^-^-(«) (2.30) 



and can be used to trade TZ for the dilaton. Similarly, one can get an explicit expression for the 
Kahler potential (|2.11|) , which to our knowledge did not appear so far in the literature. Some 
details of its derivation can be found in Appendix ^ Here we just present the final result^ 



K 



M 



7^2 1 

— K{z,z)--lmFj,^ 



\ a,b "^^lo. ^^ib / 

^ E^^. / Y ~ ^""'^"^ (?a,A -Pf ReFAs) log (1 - e"2-=-)] 



+ 



1287r4 



7a " '■lb 



^ ^n^a^^^ (7a, 7b) I T^aC I ^bC'^^, 



(2.31) 



where we abbreviated 



^^aC = y log(l-e-2-^-K)^ 



)) . (2.32) 
As we shall see, both quantities, ( |2.30|) and ( |2.31| ), appear naturally also in the context of TEA 



in the next section. Besides, note that since is symplectic invariant, the twistor lines ( 2.28 



form a vector under symplectic transformations, the contact potential is invariant, whereas the 
Kahler potential can be shown to be invariant up to a Kahler transformation (see Appendix 
^). Thus, the whole construction respects the symplectic symmetry. 



^In our normalization the Kahler potential differs from the one in p9^, nSl by a constant i?-dependent factor. 
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3. Moduli space geometry and TBA 



3.1 Relation to TBA 



As was noticed in |T3| (Appendix E of the revised version), the equations ( p.25| ), which encode 



the geometry of the instanton corrected moduh space, turn out to coincide with the equations 
of TBA. These are equations for an integrable system of particles in 1 + 1 dimensions. The 
particles are characterized by spectral densities ea{0) considered as functions of the rapidity 
parameter 6 which defines their two-dimensional momentum. TBA imposes the following non- 
linear integral equations on the spectral densities |jl2[ : 

ma(3cosh9 = eaiO) + ^ dO' cp,,{e - 9') log (l + e^^''-^''^^')) , (3.1) 

where /3 is the inverse temperature, rria are mass parameters for the particles of type a, fia are 
their chemical potentials, and (pab{0) = — i^^^^ is defined by the two-particle S*- matrix Sab{G)- 
To establish a relation between ( p.25| ) and (|3?1|), we set 

6,(e) = 27ri(S,„(ie''^"+^)-e,), (3.2) 

where ipa = argZ^^. This implies that for every charge one changes the coordinate on CP^ as 
( = 16^^^"+^ so that the BPS ray i^^ is mapped to the real axis. Then plugging these changes 
into ( p.25| ), one obtains 



oo 



W|Z,Jcosh^ = e.(^)-^5^n,,(7.,76) / dO' + log (l - e'^-Q^-^^^')) , (3.3) 

' — oo 

where ^ah = e'^^^-^-\ Comparing (§]|) with (^, one finds that they have the same form if 
one identifies 

1 + ^ ab 

I3ma = 4:'nn\Z^X /3/ia = -STTiGfe + 7ri, 0^^(6') = -- (7^, 7^) — — (3.4) 

2 e'^ - ^>ah 

whereas the additional factor n^^ in the sum is considered as a weight of the particles of type h. 
3.2 Potentials, free energy and Yang— Yang functional 



The curious relation observed in [|T3| and presented in the previous subsection can be deepen 



by considering some quantities playing an important role in TBA and comparing them with 
potentials of quaternionic geometries of the moduli spaces, the contact potential ( |2.30| ) and the 
Kahler potential ( p. 311) . 



First, the most important quantity, which is usually considered in the context of TBA, is 
the free energy of the integrable system. It is given by 



J'(/3) = ^ ^« / ^ log (1 + e'^'^--^'^^^)) 



(3.5) 
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It is straightforward to check that it coincides with the instanton part of the contact potential 
so that one has the following relation 



1 



■HP)- 



(3.6) 



On the other hand, the Kahler potential of the gauge theory turns out to be related to the 
so called Yang- Yang functional [|T7|. This is a functional which generates the action principle 
for the Bethe equations following from it by varying with respect to the spectral densities. It 
can be conveniently written as p 

a,b a 

(3.7) 

Here fai^) is the interacting part of the spectral density and Xa{0) encodes its free part together 
with the chemical potential 

ipa{0) = ea{0) -rriaP cosh 6, Xa{6) = P {fia-rria cosh 6) -ni. (3.8) 

Varying ( |3.7|) with respect to and pa and using (|3.8| ), one indeed gets the equations of TBA 
( p.l|) . The critical value of the Yang- Yang functional is then given by 




Comparison of this expression with the exact Kahler potential (|2.31|) of the gauge theory moduli 
space reveals that it reproduces two instanton symplectic invariant terms of the latter. Thus, 
one has the following relation 



where iV^^ is the inverse of N^y. = — 2ImFAE, Qab is constructed in terms of charges as 

Qab = \ NAj^p^pf + N^^ {qa,A - Re Fab) {qb,^ - Re Feh) , (3.11) 

and we used the holomorphic coordinates wa defined in ( |A.1(]| ). Although we found that the 
Yang- Yang functional does not coincide with the full instanton contribution to Km, it captures 
the most non-trivial part of the Kahler potential. In general, since Km is subject to Kahler 
transformations, the exact equality should not be expected and the appearance of the additional 
terms is not surprising. 



(3.10) 
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4. Integrable structure of instanton contributions 



4.1 The S'-matrix 



The S'-matrix corresponding to TBA 
( p.4j) . In this way, one finds 



can be easily obtained by integrating 0ab(^) from 



Sab{0) — C, 



ab 



sinh ( -{e + i{^a-^b)) 



{la rib) 



(4.1) 



where Cab is an integration constant. As was noticed in [|13|, this S-matrix is non-unitary. 
However, it is not necessarily a problem since nowadays there are many non-unitary integrable 
models. On the other hand, integrability and consistent physical interpretation require the 
S'-matrix to satisfy a set of severe conditions. In this section we show that our S'-matrix ( [4.1D 
fulfils all of them. Checking these properties, it will be important that the angles ipa associated 
to antiparticles differ from the particle angles as ipa = ipa^ tt. This relation is clear from Fig. 
Sb. 



The conditions imposed on a two-particle S-matrix include (see, for example, p8|, P9|): 



Lorentz invariance — It simply means that the S'-matrix depends on the rapidity differ- 
ence of two particles = 9a — Oh, what is clearly true in our case. 

Zamolodchikov algebra — It means that the particle creation operators must satisfy 
<l>a(0)$fe(6'') = Sab{d - d')^b{d')^aid). Applying this identity twice, one gets the following 
restriction on the S-matrix 

SabiO)SU-0) = 1. (4.2) 

In the case of unitary theories, this relation can be seen as a combination of two conditions: 
unitarity and Hermitian analyticity |2D|. As easy to see, our S-matrix satisfies only the 
combination ( [4.2|) as soon as the integration constants are chosen so that 



CabCba 



\{larib) 



(4.3) 



Crossing symmetry — It relates the scattering in the s- and t-channels and requires that 



Su{T^l-d) = Sab{e). 



(4.4) 



Again the matrix ( |4.1D fulfils this constraint if 

Cm = (-l)<^"'^''>a 



(4.5) 



• Yang-Baxter equation — It means that the order in which the particles are scattered does 
not matter and can be depicted as follows 
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ah c a be 




Sab{O)SU0 + 0')SU0') = SU0')Sac{0 + 0')Sab{0), (4.6) 

When particles do not have additional degrees of freedom, as in our case, the S matrix is 
purely diagonal and this equations is trivially satisfied. 

• Bootstrap identity — This is the most non-trivial requirement on the S'-matrix which 
relates its singularity structure to the spectrum. Namely, it demands that if Sab{0) has a 
pole^ in the physical strip, i.e., at 9 = iu'^j^ where u'^^ G (0,7r), then the spectrum should 
contain the bound state c with the mass 

ml = 1711 + 1711 + 2mam5COSM^ft, (4.7) 

appearing in the fusing process a + b c. The crossing symmetry then leads to the 
existence of other two fusing processes, 6 + c — )■ a and c + a — > 6, with the fusing angles 
satisfying 

u:, + ul + ul = 27r. (4.8) 

But the most important condition is that it does not matter whether an additional par- 
ticle, say d, scatters with the bound state c or consequently with the two particles a, b 




S,ai0 - ml)S,,ie + ml) = S,,ie), = vr - ul,. (4.9) 

In our case the poles of the S'-matrix ( [4.1|) correspond to 

U:, = A-^a. (4.10) 



^Usually it is taken in the strong form of a pole of order one, but experience in the field of S matrix shows 
that this principle has to be usually extended to poles of order greater than one Q . 
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The condition that the pole is on the physical strip is equivalent to sinw^^ > 0. The mass 
formula (|4.7| ) together with its expression ( p.4|) in terms of the central charge yields 



/3m5 = 47^7^|Z^,+^J, (4.11) 

so that we deduce that the bound state c has the charge 7a + 76 consistently with the 
physical interpretation. Moreover, since (76,7c) = {ic^la) = {la, lb), the S-matrix ele- 
ments Sbc and Sea have poles of the same degree as in Sab- By simple manipulations one 
finds that the corresponding fusing angles® u^^ = "V^c ~ i^b, u^ca = "^a — '^/'c + 27r also belong 
to the physical strip due to 

sin< = — sinM^„ = — sin<fe > (4.12) 
mb rric 

and satisfy the constraint ( |4.8| ). This means that the other two fusing processes obtained 
by crossing also exist. Finally, it is straightforward to check that the bootstrap identity 
I) does hold provided 

CdaCdb = Cdc- (4-13) 



The three conditions ( [4.3|) , (|4.5| ) and (^4.13 ) on the integration constants fix them to be 



Cab = a^feC^^'^") where aab is a Z2-valued function on the square of the charge lattice 
satisfying the same conditions ( [4.3|) , ( |4.5D and ([4.13|) . Given a polarization into electric 
and magnetic charges, it can be chosen, for example, as aab = (—lY'^'^Pb, It reminds a 
lot the quadratic refinement introduced in [|T3] and mentioned in footnote ^. 

Thus, our S'-matrix satisfies all necessary requirements of an integrable model in 1 + 1 
dimensions. 

4.2 Y-system 

Y-system pi], is a system of functional algebraic relations on the exponentials of the spectral 



densities which, although equivalent to the integral equations of TBA, play an important role 
in the analysis of integrable models. In particular, it may be viewed as an intermediate step 
between TBA and the transfer matrix approach. Therefore, it would be nice to have such 
relations at our disposal. 

Let us introduce the Y-functions as^ 

Ya{e) = -e^i^^-'-io) , (4.14) 
Using ipa = i'a + T^ and 9^ = —6a, it is easy to check that ( p.3| ) leads to the following relations 

Y. (O + f ) n (« - = n [l + ' (f + - (4.15) 

This is our Y-system. Of course, for a general configuration of charges, it is extremely com- 
plicated. We would like to point out however a few unusual features of this Y-system (|4.15|) 
comparing to the standard Y-systems appearing in the literature on integrable models: 

^These definitions are valid for i/^a < i/Jc- Otherwise the shift 2tt will appear in one of the other fusing angles. 
^Note the additional minus sign in the definition of the Y-function. It is introduced to ensure that = Y" 
and hints that the bosonic version of TBA might be more relevant in this context (see the end of section H). 
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First, on the l.h.s. of ( 4.15|) one multiplies functions associated with a particle and its 



antiparticle, whereas usually one has only one function. This feature is related to the 
absence of the parity symmetry in our case and to the unusual reality conditions. Notice 
that the reality conditions ( p.27| ) in terms of the spectral densities and the Y-functions 
read 

eaiO) = e-ai-e), Yaie) = Yai-e). (4.16) 

As a result, the Y-functions are not necessarily real on the real axis of the spectral 
parameter. Neither is the combination on the l.h.s. of ( |4.15|) . 



Second, on the r.h.s. the Y-functions are all evaluated at different points, whereas usually 
their arguments do not contain any shifts. Moreover, usually the fusing angles and the 
only shifts appearing in functional relations are rational multiples of vr. Here they are 
completely arbitrary and vary continuously with the moduli 

Third, usually the power of each element of the product on the r.h.s is related to the 
incidence matrix of a graph which structure is severely constrained by the periodicity 
of the Y-system |S^, This incidence matrix is equivalently described by the ma- 
trix Nab = 4>ab{d)d9. In our case it is just not defined because the kernel is not 
integrable. However, this describes only a particular class of the Y-systems which are 
known nowadays. It is known that the general mathematical structure which lies behind 
Y-systems is related to cluster algebras [^, |3^. Such general Y-systems are defined 
by skew-symmetrizable matrices having to satisfy some mutation identities. Our matrix 
'^7i,(7a,7b) is also skew-symmetrizable and originates in some mathematical structures 



also having relations with cluster transformations [^. Therefore, it would be quite in- 
teresting to understand what is the precise connection between cluster algebras and the 
Y-system appearing in our case. 

4.3 Remarks on the strong coupling limit 

Usually, TEA is a very effective tool to get the conformal or high temperature limit of the inte- 
grable model. This is a limit where the parameter /3 ( ^.Ij ) goes to zero. Given the identification 
( |3.4]) , this is equivalent to vanishing of the parameter TZ, which means either the small radius 
limit in gauge theory, where it becomes effectively three-dimensional, or the strong coupling 
limit in string theory, where the ten-dimensional string coupling becomes large. 

The standard derivation tells us that in this limit the Y-functions become constant and 



real for a wide range of the rapidity parameter 6 ||T^ . The values of these functions can then 
be easily found from the Y-system which reduces to a system of algebraic equations. Finally, 
there is a nice formula 

7r ^ — ' 

a 

which gives the free energy in terms of these constant values, where L is the Rogers dilogarithm 
defined by C{x) = Li2(x) + | log(a;) log(l — x). Thus, if the same story was valid in our case, we 

^°It is interesting that there is a case where many of these comphcations disappear. We present it in appendix 
^ This is probably one of the simplest possible examples of the TBA systems of the type considered in this 
paper. 
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could hope to find the strong couphng hmit of, for example, the contact potential in an easy 
way. 

However, the TBA (|3.3| ) has several features which distinguish it from the usual integrable 
models and make the story much more complicated. In particular, one has: 

• more complicated reality conditions ( [4.16|) than the ones appearing usually, 

• the kernel 4>ab{0) (|3.4| ) is not decaying at infinity, 

• the system is supplied by arbitrary imaginary chemical potentials fia- 

On the other hand, the standard derivation of the free energy in the conformal limit relies on 
the absence of these features. In fact, it can be extended to include some of them. For example, 
p6[ gave a generalization for non-vanishing chemical potentials leading to the same formula 
( ^.17| ) with a simple modification of the Rodgers function C: in logya in the second term one 



should subtract the contribution of the chemical potential, i.e., it can be replaced by — e^. But 
altogether the above features give rise to the appearance of new phenomena. Most importantly 
is that the Y-functions are not constant anymore in this limit. This can be seen analytically 
and has been also verified by a numerical analysis of the simplest example proposed in appendix 
p. As a result, the derivation of the free energy must be seriously reconsidered and we leave 
the detailed investigation of this issue for a future work. 

There is also another problem to handle in order to find the strong coupling limit of the 
free energy, which comes from the fact that the lattice of charges is infinite. Indeed, the formula 
( ^.17| ) gives the free energy as a sum over all particles in the spectrum so that we will have to 
sum over all charges. This sum is usually divergent and requires a certain resummation P? 



But the resummation should be performed before the limit since these two procedures are not 
commuting. 

This issue can easily be exemplified if one considers the contributions of only D(-l)- 
instantons to the hyper mult iplet moduli space of string theory. In this case TBA can be 
solved exactly since the scalar product of two charges is always vanishing leaving us with 

eq(^) =47^7^|g|coshe. (4.18) 

Then plugging this result into (|3.5| ) and extracting the limit or using ( [4.17|) with Yg = e~'^'"^'^'^° , 
one finds (with Ug = xx) 

HO) = -^J2^i2 {e-'''''"^') (4.19) 

which is clearly divergent. On the other hand, the same limit can be found by a Poisson 
resummation of the initial expression obtained by taking the integral in ( |3.5|) explicitly in 



terms of Bessel functions |38, 15 1. This leads to p3 



•^(/3) = -ff C(3)r' + 0(/3), (4.20) 
where we identified f3 with TZ. The leading contribution has a different scaling in /3 than the 



one of the usual result (|4.17|) and reflects the divergence of the latter. 

Thus, the study of the strong coupling limit is supplied here with two problems: first, 
our TBA is much more general than one usually considers and, second, the infinite spectrum 
requires a resummation technique to be applied. 
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5. Discussion 



In this paper we demonstrated that the relation noticed in |13| between the equations describing 
the non-perturbative geometry of moduh spaces and the equations of TBA goes beyond the 
formal analogy. To this end, we provided an identification of the physically relevant quantities 
at the two sides of the correspondence, showed that the S-matrix underlying this TBA fulfills 
all usual constraints imposed by integrability, and derived the associated Y-system. 

Of course, these results are only a first small step towards a deeper understanding of this 
relation. However, already at this point we seem to open new interesting connections with other 
developments. First, quite similar TBA equations and Y-systems to those considered in this 



paper appear in the context of minimal area surfaces in AdS^ |10]. In fact, the TBA for the 
full AdS^ problem found in |^ possesses all features listed in section [4.3| which complicate the 
evaluation of the conformal limit. Nevertheless, it has been successfully computed in [Q. This 
was possible due to an additional Z4 symmetry, which considerably simplifies the corresponding 
Y-systems, and additional restrictions on the chemical potentials. The TBA with generic 
chemical potential seems to experience the same phenomena which were mentioned in section 
O] and thus requires much more care. 

Second, the Yang- Yang functional, providing for us the Kahler potential on the gauge 
theory moduli space, plays also a prominent role in the relation of BPS vacua of certain N = 2 
gauge theories to Bethe states of integrable models p. In that correspondence it has been 
identified with the twisted effective superpotential of the low energy effective theory. Since the 
context of the two stories is very similar, it is tempting to assume that the Kahler potential 
and the twisted effective superpotential are also related. 

Besides, a somewhat intriguing observation is that the Y-system presented in appendix ^ 
leads to the MacMahon function, the generating function of 3d partitions. This hints that there 
might be a relation with the beautiful duality between melting crystals and topological strings 



A set of interesting questions arises if one tries to draw physical consequences from the 
fact that the S'-matrix ( [4.1| ) is consistent with the bootstrap. Usually, the latter can be used 
to generate the complete spectrum out of some "elementary" particles. Can it be used in the 
same way, for example, for gauge theory? Which singularities of the S'-matrix should be taken 
into account, all or only those which belong to the physical strip, as in the main text? Do zeros 
of the S'-matrix play some role? 

The last two questions become especially relevant, if one remarks that there are actually 
two versions of TBA, for fermionic and bosonic particles. Here we used the fermionic version. 
The bosonic equations differ only by two signs: one should flip the sign in front of the integral 
in and the sign in the logarithm inside the integral. This does not change much. The 
second sign leads to disappearance of the shift by vri in the chemical potential fia ( p.4|) . But the 
first sign gives the S'-matrix which is the inverse of ([4.1|) . As a result, the zeros and poles are 
exchanged between each other. What version of TBA is relevant for our problem? Is it related 
to the signs of the invariants n-y or to the quadratic refinement of |]13|? Is it important at all? 
These are just few questions which arise naturally. 
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A. The Kahler potential 



In this appendix we derive the exact instanton corrected Kahler potential (p.31| ) and after that 
discuss its transformation under symplectic symmetry. 

A.l Evaluation 

The starting point for evaluation of the Kahler potential is the representation ( ^.101 ), which in 
our case takes the following form 



K 



M 



a 



Li2 (f 



O + ^^-^-APa (log (1-e 



+27rip^ Ul^-B^ log(l-e-2-=-) 



(A.l) 



where = qa,AV\^] "Pa/^lv'' twistor lines can be read off from ( p.28|) , ( |2.16|) , and 

the coordinates A^, B\ are related to the physical fields as follows |1^ 



A' 



(A.2) 



The expressions for r^^j will not be needed explicitly. For our purposes it is enough to know 
that 

<±](C) = r/f,](C) + 0(C±2). (A.3) 

Plugging all definitions into ( |A.1| ), evaluating the first two integrals by residues and taking into 
account the quasi-homogeneity property (|2.14| ), it is straightforward to obtain that the Kahler 
potential is given by the first two lines in ( p. 311) plus the following contribution 



skr^T.^'^^Pa [ T^aC n,„ga,Alog(l-e-2-=-)-^^n,,g,,A /" V,C 



c-c 



(A.4) 



with the measure VaC defined in (|2.32|) . First, let us concentrate on the second term for h ^ a. 
The double integral is antisymmetric in a, h so that it reproduces the remaining contribution to 
Kji^ from the third line of ( |2.31|) . Second, the term with h = a seems to be singular. However, 
it is easy to realize that it can be represented as a difference of two double integrals such that 
the first integral goes over a contour which is either to the left or to the right of the second. 
As a result, it gives just the residue at C' = C which exactly cancels the first term in (|A.4|). 
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A. 2 Symplectic invariance 



The Kahler potential on the moduh space is expected to respect the symplectic invariance 
of the gauge theory. However, the potential (|2.31|) is clearly not invariant under symplectic 
transformations. This is true already at the tree level. This phenomenon was explained in 



40| where it was shown that is actually invariant up to a Kahler transformation so that 
the moduli space metric does not change. Here we would like to generalize this result to the 
non-perturbative level. 

To this end, let us review how the symplectic transformations affect various quantities 



411 These transformations are represented by 2d x 2d matrices acting on the symplectic 
vectors as 



whose blocks satisfy 



A'B - B^A = J^C - C^A = B^C - C^B = B^V - V'^B = 0, 

A^V - C^B = 1. 



(A.6) 



It is convenient to introduce the holomorphic matrices 

5^(2) = + i3^®FeE(^), Z^^{z) = [S'\z)]\B''''. (A.7) 

The nice feature of the matrix Z is that it is symmetric whereas S encodes transformation 
properties of various quantities such as 

i^AE ^ (I?A=i^=e + CAe)[5-i]®2, 

(A.8) 

S' 



Ca - F^^e ^ (Ce - i^eEC"") [S~'f 
QA - Faep"" ^ {qe - Fe^p^) [S-^f 



Using these results, one can show that the variation of the Kahler potential ( p.31| ) under a finite 
symplectic transformation is given by 

AKm = -^Z^^(^)w;a«;e +'-Z''^{z)wAWj,, (A.9) 



where 



together with z^ = — iz^[^]|(=o provide the complex coordinates on the moduli space Ai. This 
demonstrates that the change of Km is described by a Kahler transformation in agreement 
with symplectic invariance. 
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B. Example: rigid Calabi-Yau 



Let us consider a particular simple case of a compactification on a rigid Calabi-Yau X with 
the vanishing Hodge number h2^i{X). Then the hypermultiplet sector consists only from one 
hypermultiplet known as the universal hypermultiplet. As a result, the lattice of charges is two- 
dimensional, 7 = {q,p)- In addition, the holomorphic prepotential is simply F{X) = — 
where r is a fixed complex coefficient determined by the holomorphic 3-form of the Calabi-Yau 
121. As a result, the central charge function becomes 



= q + rp. 



(B.l] 



Let us restrict ourselves only to the two sets of charges: pure "electric" (g, 0) and pure 
"magnetic" (0,p). In other words, we ignore all possible "dyons" with both electric and mag- 
netic charges non- vanishing. From ( |B.1| ) it is clear that in this case there are only four sets of 
angles in the game {q,p > 0)^^ 



^(q,0) = 0, '0(O,p) = A 



vr. 



(B.2) 



where ipr = argr. Since Yg^ = Y^q and Fq.p = ^o'l, in fact one has only two unknown functions. 
As a result, in this sector the TBA equations read as follows 



47^7^ cosh^ + 27riC 



no,pP 



p>0 



de' 



logYi,oW 

e' 



^"-,\og[i-Yi,{e'))-'-, 



AttU cosh 9 - 27riC = - log Fq,! (^) 



477 ^ 



q>0 



+ r 



log (1 - Yi,{e')) 



^ log (1 - Yl,{-Q')) 



-7 log (1 - n'!o(-^')) 



(B.3) 



where we defined r' = e^"^^ . The corresponding Y-system ( [4.15|) takes in this case the following 
form 



^1,0 



TTl 

y 



TTl 

vri 
2" 



n i-^cTi (^+i(f-v^o) 



p>0 



n[i-no (^-i(f-^o) 

g>0 



qUqfi 



(B.4) 



Especially simple the above equations become when the parameter r is pure imaginary, 
i.e., Vr = 7r/2. Then, all angles ( [B.2| ) are multiples of 7r/2 and all fusing angles are multiples 
of TT and moduli independent, as in the usual integrable models! Moreover, in this case the 



^^One could ask whether the restriction to only electric and magnetic charges is consistent with the bootstrap, 
which as we know leads to the bound states of charges 7a + 7h. Remarkably, for Imr > 0, all fusing angles 
corresponding to the poles of the S'-matrix (4.1) with '0a from (B.2) turn out to be outside of the physical strip. 
Therefore, if such poles are not required to satisfy the bootstrap identities, our restriction is consistent. 
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Y-system loses one of its unusual features. Namely, in (|B.4| ), as in the standard Y-systems, all 
functions on the r.h.s. appear without any shifts in the arguments. 

To analyze the resulting Y-system, it is useful to take into account that Ug^ = xx 



One may assume that no,p is also p- independent. Then the r.h.s. of ( B.4 ) is described by the 
MacMahon function^^ 



n=l 



(B.5) 



This is the generating function of plane partitions which describes also the large volume limit 

This fact played the 
Remarkably, 



of the topological string partition function as Z{gs) = [S{e 9')'\^x/2 



prominent role in the duality between topological strings and melting crystals [18 



in our example not only this function appears, but even the power is also given by the Euler 
characteristic of the Calabi-Yau. This suggests that there might be a deep interplay between 
our story on one side and random partitions and topological strings on the other. 

To exemplify some properties of the TBA systems introduced in this paper, it might be 
useful to consider a further truncation of this example, where one drops the sum over the 
infinite set of charges and considers, for example, only charges with q = ±1 or p = ±1. Such 
a truncated system is particularly suitable for numerical analysis, but still possesses most of 
the non-trivial features of the full problem such as non-decaying kernels, arbitrary chemical 
potentials, non-trivial reality conditions, etc. 
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